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We describe the resonance fluorescence spectrum of an atomic three-level system where two of 
the states are coupled by a single monochromatic laser field. The influence of the third energy level, 
which interacts with the two laser-coupled states only via radiative decays, is studied in detail. For a 
suitable choice of parameters, this system gives rise to a very narrow structure at the laser frequency 
in the fluorescence spectrum which is not present in the spectrum of a two-level atom. We find those 
parameter ranges by a numerical analysis and use the results to derive analytical expressions for 
the additional narrow peak. We also derive an exact expression for the peak intensity under the 
assumption that a random telegraph model is applicable to the system. This model and a simple 
spring model are then used to describe the physical origins of the additional peak. Using these 
results, we explain the connection between our system, a three-level system in V-configuration 
where both transitions are laser driven, and a related experiment which was recently reported. 

PACS numbers: 42.50.Lc, 42.50.Ct, 42.50.Hz 



I. INTRODUCTION 



There have been numerous contributions to study the 
various aspects of the resonance fluorescence of laser- 
driven few-level atomic systems. One topic of particular 
interest is the appearance of narrow lines in the fluores- 
cence spectra of laser-driven three-level systems. In these 
systems, the additional third atomic state together with 
its coupling to the other two states allows for substantial 
changes in the resonance fluorescence spectrum, as com- 
pared to the spectrum of a two-level atom In Q the 
appearance of sub- natural line widths was predicted for a 
three-level system in V-configuration, a prediction which 
was verified experimentally |3j. Also in three-level sys- 
tems in A - configuration the possibility of line narrowing 
was found Q . In j^] , the spontaneous decay between two 
atomic levels was modified by an additional atomic level 
coherently coupled to one of the other states. One of 
the results was a possible occurrence of narrow lines in 
the fluorescence spectrum. Q found narrow lines due to 
quantum interference in an atomic system with two close 
energy levels which couple to the same vacuum modes. 
Line narrowing can also be found in systems with more 
than three energy levels. reports spectral narrowing 
controllable by the phase difference between two driv- 
ing laser fields in a four level system. In principle it is 
also possible to have unbounded line narrowing in the 
emission spectrum with essentially complete transfer of 
the corresponding intensity into the narrow feature |^J. 
Somewhat similar, also the probe absorption of an atomic 
few-level system may be altered favourably under the in- 
fluence of driving laser fields (See e.g. M). 
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A different mechanism for the appearance of narrow 
lines is electron shelving [|l0| [ll]. In 1 12, this phe- 
nomenon was discussed in a three-level atomic system 
driven by two laser fields in a V - configuration. One of 
the fields is assumed to be much weaker than the other. 
This allows for a shelving of the atomic population into 
the weakly coupled atomic level. Because of the long life- 
times of such shelved states, the fluorescence spectrum 
may exhibit narrow peaks. Compared to the usual two- 
level fluorescence spectrum, the spectrum of the atomic 
system discussed in |l2|, |l3| contains an additional very 
narrow peak centered at the laser frequency for a suitable 
choice of parameters. The peak is controllable to some 
extent by altering the strength of the weak coupling laser. 

Recently there has been an experiment testing the 
ansatz of electron shelving In the experiment a 

three-level atomic system different from the one de- 
scribed above was used, in which only one transition is 
driven by a laser field. The third atomic level was cou- 
pled to the driven levels by radiative decays as shown in 
Fig. |l|. Experimentally, this was realized by a setup as 
in Fig. ^|. In the experiment, the additional narrow peak 
due to an electron shelving was found. But although the 
atomic system used was different, the experimental re- 
sults were reported to be consistent with the theoretical 
results in ft^ |. 

In this paper, we provide the theoretical background 
for the atomic system used in the experiment mentioned 
above. We set up the equations of motion and calcu- 
late the resonance fluorescence spectrum. By a numer- 
ical analysis, we find suitable parameter ranges for the 
exhibition of the additional narrow peak, which may be 
useful for further experiments. We use these results to 
derive analytical expressions for the additional peak in 
an approximation which is similar to the secular limit for 
high laser intensities. Also we derive an exact analyti- 
cal expression for the peak intensity valid for any choice 
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FIG. 1: The system of interest in bare state representation. FIG. 2: Possible experimental realization of the atomic sys- 
fi is the Rabi frequency of the driving laser field, and 27, 272 tem of interest. R represents a pumping from j2) to |4). 
and 273 are the total decay widths. 



of parameters under the assumption that a random tele- 
graph model is applicable to describe the physical origin 
of the narrow peak. Then we use this random telegraph 
model to clarify the connection between the experiment 
and the two different atomic systems in a discussion of 
the physical origins of the additional narrow peak. Fi- 
nally, we present a simple spring model to give an intu- 
itive picture of the atomic dynamics. This model is able 
to reproduce the different components found in the res- 
onance fluorescence spectrum and allows to easily derive 
approximate expressions for their widths. 

The paper is structured as follows. In the second sec- 
tion we derive expressions for the resonance fluorescence 
spectrum of the given atomic system. In section three, we 
use the derived equations for a numerical analysis of the 
resonance fluorescence spectra. The results of this analy- 
sis motivate the approximations made in the first part of 
the following analytical treatment of the problem. In the 
second part of the analytical section, we use the telegraph 
model to derive expressions for more general parameter 
ranges. In section four we give two different physical in- 
terpretations for the occurrence of the additional narrow 
feature in the fluorescence spectra, the random telegraph 
model and a spring model. The first model is used to 
clarify the connection between our system of interest, the 
system where both transitions are driven coherently and 
the experiment. The second model is the spring model 
for an intuitive picture of the atomic dynamics. Section 
five will discuss and summarize the results. 



II. PRELIMINARY CONSIDERATIONS 

In this section, the atomic system studied in this paper 
is described. We set up the equations of motion and use 
these to derive expressions for the resonance fluorescence 
spectrum. 

A. System of interest 

The atomic system considered in this paper is shown in 
Fig. [|. hwi (i — 1,2,3) are the respective energies of the 
three levels, and 2j, 272 and 273 are the decay constants 



(Einstein coefficients). The transition 1 <-> 3 with the 
frequency luq — (073 — u>i) is driven by a laser with the 
frequency u>l and the Rabi frequency Q. A = (lql — u)q) 
is the detuning. In the analytical calculations, we take 

72,73 < 7 (1) 

as a general assumption as this will turn out to be the 
most interesting parameter range. We neglect all terms 
of order higher than 74/7 (i = 2, 3) in the analytical cal- 
culations if not noted otherwise. Thus the second atomic 
level is only weakly coupled to the other levels. In an 
experiment this is either possible by using a quadrupole 
transition for 3 «-> 2 or 1 <-> 2 or by using dipole transi- 
tions including a fourth ancilla level as shown in Fig. |^. 
In this figure, R is a pumping mechanism which transfers 
atomic population from |2) to |4). 

B. Dynamical equations 

The Hamiltonian H for the system is 

H = Hp + Hi 
H F = a(o; 1 -A)|l)(l|+^|2)(2| + Rw3|3)(3| 

Hj = fiA|l)<l|-^(|l)<3| e ^*+ h.c.) 

The Hamiltonian V in the interaction picture with re- 
spect to Hf is 

K = fiA|l)<l|-^(|3)(l| + |l)(3|) (2) 

The equations of motion for the corresponding density 
matrix p in the interaction picture are 



P33 


= -2(7 4 


- 73)^33 - 


^ y(P3i - P13) (3) 


Pll 


= 27^33 - 


- 272P22 


-y(P3i~Pi3) (4) 


hi 


= y(P33 


- Pn) - 


- (7 + 73 ~ iA)p 3 i (5) 


> = 1 


—psa—pii 


and pij 


= ( Pji )*iovi,je {1,2,3}. 



P32 and P21 relax with the rate (7 + 73 +72) and 72 
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respectively. Introducing a = 72/73, the steady state 
solution p ss to this system becomes 

P33 = Jj aii 

pit = ^a(4A 2 + 47 2 + 8 7 73 + 473 2 + ^ 2 ) 
P22 = ^ 2 

Psi = ~^2ail(A + i( 1 + l3 )) 

P32 = P21 = (6) 

with the abbreviation 

N = fl 2 + 2a(2A 2 + 2 7 2 + 4773 + 2 7 f + fl 2 ) 

Upon inserting a = 72/73 and setting 73 to zero in the 
above expressions, the result simplifies to the usual ex- 
pressions for a laser-driven two-level system. Throughout 
the paper, a superscript "2-level" denotes values for this 
two-level system, whereas a superscript of "3-level" or 
no superscript means the corresponding values for the 
three-level system discussed here. 



C. Observables of interest 

The equations of motion Eqn. (§)-(||) can be written 



as 



P = (P33,Pll,P31,Pl3) J 



a = (<Ji,(J2,cr 3 ,(T i ) T 

= (|3)(3|,|l)(l|,|3)(l|,|l)(3|f 

with time independent time evolution matrix B and a 
constant vector / given by 



( -2 (7 + 73) 
2 (7 - 72) -2 72 



V 



-m -iA-7-73 







iA-7-73 / 



and/= (0,2 72 ,0,0) T . 

The normalized resonance fluorescence spectrum S(lo) 
is then given by the real part of the Fourier transforma- 
tion of the two-time correlation function {a^^a^t')) as 



— Re / < 

TTPi Jo 



-^-^{a^Oiit')) dtdt' (7) 



The normalization is chosen such that the integral of 
S(lu) over all frequencies lj yields unity. Thus the nor- 
malized expressions for the different contributions to the 



total resonance fluorescence spectrum are equal to the 
relative intensities of the spectral components. 

Following [plU , we split the expression for a into its 
average value and the deviation from the average 
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for i £ {1, 2, 3, 4}. By inserting this, the correlation func- 
tion can be separated into two terms corresponding to its 
average motion and its fluctuations around the average: 

(<rs(t)a*(t')) - (a 3 (t))(a 4 (t')) + (Sa 3 (t)Sa 4 (t')) (8) 

The first term corresponding to the average motion yields 
the coherent or elastic part S co h{uj) of the normalized 
spectrum 
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Scoh^) = ^-6(u-U L ) 

P33 



Thus we have 



'coh 



Pi 



(9) 



(10) 



as the relative intensity of the elastic peak. 

The second term represents the incoherent or inelastic 
part of the spectrum. According to the quantum regres- 
sion theorem we have for t' > t and i £ {1, 2, 3, 4} 

^{a(t')a i (t))=B{a(t')a i (t))+I{a i (t)) (11) 

This together with Eqn. ([?]) and (||) yields for the total 
normalized incoherent spectrum Si nc {uj) consisting of the 
Mollow spectrum and the additional narrow peak 
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Re 



KP33 
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(12) 



where the brackets with the subindex mean that the 
fourth component has to be taken and "Re" indicates 
that the real part of the following expression has to be 
taken. Also we have 

R = (fo^O)^ (0)) ss 

= (<?(o)<r 4 (o)r - <<?(o)) s >4(o)r 

where the index "ss" denotes the steady state value. 



III. CALCULATION OF THE SPECTRA 

In this section we calculate the resonance fluorescence 
spectra emitted by the given atomic system. First the 
equations derived in the previous chapter are used for a 
numerical analysis without approximations. In the sec- 
ond part, we derive analytical expressions for one of the 
spectral components using the results of the numerical 
analysis. 
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FIG. 3: Sample spectrum. Parameters used are 7 = 1, 
Q = 6, 73 = 0.005, a = 0.3, A = 0. The solid line shows the 
full normalized inelastic spectrum. The dotted line shows the 
normalized Mollow spectru m extracted from the full spectrum 
as shown in section 
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FIG. 4: Closeup on the additional narrow peak. Plots and 
parameters are as in Fig. |^. Note the very different axes 
scales. 



A. Numerical analysis 

The sum of the elastic part in Eq. (||) and the inelastic 
part in Eq. ( |l2| ) give the full normalized resonance flu- 
orescence spectrum without approximations other than 
those adopted in the derivation of the equations of mo- 
tion. In this section we use this expression for a numerical 
analysis of the spectra emitted by the system of interest. 
In all plots, the elastic component is omitted. Thus nar- 
row peaks centered at the laser frequency are always new 
structures due to the addition of the third atomic level. 

A sample spectrum is shown in Fig. [| The solid line 
shows the full normalized inelastic spectrum. The dot- 
ted line shows the normalized Mollow spectrum extracted 
from the full spectrum. As elaborated later, this Mollow 
spectrum is the spectrum one would expect from a driven 
two- level system jlj . The additional narrow peak at ojl 
can easily be seen. A closeup on this peak with very 
different scales on both axes is shown in Fig. |J. The 
additional peak can be distinguished from the coherent 
peak due to its finite width. 

Since the elastic and the additional narrow peak are 
both centered at the laser frequency, experimental obser- 
vation of the narrow peak is a challenging matter. In re- 
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FIG. 5: Normalized inelastic spectrum for different values of 
73. a — 72/73 is kept fixed at a = 0.3. The other parameters 
are 7 = 1, ft = 8, A = 0. (a) 73 = 0.05; (b) 73 = 0.01; (c) 
73 = 0.005; (d) 73 = 0.001. The z-axis has arbitrary units 
and a pedestal due to the Mollow peak was removed. 
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FIG. 6: Normalized inelastic spectrum for different values 
of O. The other parameters are 7 = 1, 73 = 0.01, a = 0.3, 
A = 0. (a) Q = 10; (b) Q = 8; (c) Q. = 6; (d) ft = 4. The 
z-axis has arbitrary units and a pedestal due to the Mollow 
peak was removed. 



cent experiments, optical heterodyne detection was used 
to measure fluorescence spectra with sufficiently high 
resolution^ 
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Fig. H shows spectra obtained for different values of 
73. 72 is adjusted according to the fixed ratio a = 72/73, 
and also the other parameters are kept fixed. Effectively, 
the figure shows the influence of the coupling strength 
of the additional atomic level to the laser-driven ones 
on the additional peak. Only a very narrow frequency 
range around the laser frequency is plotted. The z-axis 
has arbitrary units, and the contribution from the Mollow 
peak centered at the laser frequency has been subtracted. 
It can be seen that for large 73, i.e. strong coupling 
of the additional atomic state to the driven ones, the 
additional narrow peak becomes broader and seems to 
vanish. For decreasing 73, the peak amplitude increases 
while its width decreases. As long as Eq. ([[]) is satisfied, 
the relative peak intensity is unaffected by the value of 73 
if a and the other parameters are kept fixed. Therefore 
in experiments it is possible to choose 73 and 72 as small 
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FIG. 7: Absolute intensity of the narrow peak component 
of the spectrum against fi and A. Parameters used in the 
plot are y = 1, 73 = 0.005, a = 0.3. The solid lines show the 
detuning calculated using the condition in Eq. . 
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FIG. 8: Sample spectrum with A = A max . The other pa- 
rameters used are as in Fig. |§| The solid line shows the 
full normalized inelastic spectrum. The dotted line shows the 
sum of the extracted normalized Mollow spectrum and the 
approximation for the narrow peak in Eq. (EOj) . 



as possible compared to 7 in order to have a narrow peak 
with large amplitude. This justifies our assumption Eq. 
(|l|) for the analytical part later. 

A comparison of spectra for different Rabi frequencies 
Q, is shown in Fig. pi Again, only a narrow frequency 
range around the laser frequency is shown, and the z-axis 
has arbitrary units with a removed Mollow contribution 
as in the previous plot. The narrow peak at van- 
ishes for increasing f2, just as the elastic component not 
shown in the picture does. This shows that both the elas- 
tic and the additional narrow component are non-secular 
features. If the Rabi-frequency is chosen too large, their 
contribution to the complete spectrum will be small. 

It is however possible to choose an appropriate detun- 
ing A to have large intensity in the narrow peak even 
for very large Rabi frequencies. This can be seen from 
Fig. [7] which shows the absolute intensity of the narrow 
peak plotted against the Rabi frequency and the detun- 
ing. For any fixed value of A, the peak intensity vanishes 
for increasing Rabi frequency f2. However by appropri- 
ately changing A with Q, it is possible to move along the 
saddle of maximum intensity in Fig. 7. In Fig. @ and 
Fig. [| compared to Fig. || and Fig. 4 only the detun- 
ing has been changed to a suitable value A max estimated 
from Fig. [7]. Nevertheless the peak amplitude has in- 
creased by more than one order of magnitude, while the 
width at half height has only decreased by a factor of 
approximately two thirds. Thus the peak intensity with 
the appropriately chosen detuning is more than 7 times 
larger than in the case with zero detuning. The exact val- 
ues for the amplitudes A, the widths T and the intensities 
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/ of the narrow peak in the plots are: 
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Closeup on Fig. fel Note the very different axes 



This is a clear indication for the fact that for zero de- 



tuning, the Rabi frequency chosen for the plots is already 
large enough to account for a secular suppression of the 
peak intensity. 

Later we will derive an analytical expression for a suit- 
able detuning A max . In the limit fl 7, 7 ^> 73,72 
and A = A max (see Eq. (|l|)), modifying SI, 7 or 73 
and 72 with constant a by one order of magnitude typ- 
ically changes the absolute intensity of the peak by less 
than one percent. The absolute intensity of the addi- 
tional peak only depends on a = 72/73- It reaches its 
maximum value for a = (2v / 3) _1 - For this value of a, the 
relative intensity of the additional peak is approximately 
24%. Therefore this limit is suitable for further exper- 
iments and might even be used to measure the relation 
between 73 and 72. 

The dependence of the narrow peak on a is shown in 
Fig. |Io|. In all the curves, 73 is kept fixed, and 72 is 
adjusted according to 72 = 073. Thus this figure shows 
the dependence of the narrow peak on the relation of the 
couplings of the additional atomic level to the upper and 
to the lower laser-driven level. The dotted line shows the 
Mollow spectrum of the corresponding two-level system 
as a reference. For increasing a, the peak width increases 
and the peak height decreases. In contrast to the in- 
dependence on 73, the peak intensity does depend on a. 
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FIG. 10: Normalized inelastic spectrum for different values of 
a. The other parameters are 7 = 1, 73 = 0.01, Q — 8, A = 0. 
(a) a = 0.01; (b) a = 0.1; (c) a = 0.3; (d) a = 1. The dotted 
line shows the corresponding two-level Mollow spectrum. 
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FIG. 11: Normalized inelastic spectrum for different values 
of a. The parameters are chosen as in Fig. [l(j. The dotted 
line shows the two-level Mollow spectrum. The four solid lines 
show spectra for the system discussed here. Apart from the 
narrow peak, the different curves are almost on top of each 
other. 



Therefore it is not possible to decrease the peak width by 
choosing very low values for a, as this would be equiva- 
lent to very low values for 72 . A low value for 72 however 
traps the population in the additional atomic state and 
thus yields a low total fluorescence intensity. The value 
a = 0.3 as chosen in most of our plots is close to the value 
which yields a maximum absolute peak intensity of the 
narrow peak. 

The Mollow contributions to the inelastic spectrum are 
almost unaffected by changing a as shown in Fig. [ll]. 
While Fig. |l0| only shows a very narrow frequency region 
around the laser frequency, this figure depicts the full 
inelastic spectrum for the same sets of parameters as in 
Fig. [l(]. The dotted spectrum for the two-level system 
and the four three-level spectra only differ in the small 
frequency region shown in Fig. [l^, where the additional 
narrow peak dominates. A similar plot shows that the 
Mollow spectrum also does not depend on the value of 
73 and 72 as long as Eq. (||) is fulfilled. 



B. Analytical considerations 



In this subsection, we derive analytical expressions for 
the spectral contributions by the additional narrow peak 
using two different methods. The results of the numerical 
analysis encourage us to focus on a restricted domain of 
parameters. Thus in the first method, we assume Eq. (^) 
and choose the Rabi frequency to be large as compared 
to the spontaneous emission widths (secular limit). In 
this parameter range, approximations to simplify the re- 
sults are possible. In the second method, we calculate 
an exact expression for the intensity of the peak under 
the assumption that a random telegraph model can be 
applied to describe the origin of the narrow peak. This 
assumption will be justified in section IV A. 



1. Secular approach 

In order to obtain an analytical expression for the addi- 
tional peak, a Lorentzian shape is assumed, because the 
usual Mollow peaks are also approximately Lorentzian 
in the parameter ranges where they are well separated. 
In the Mollow case, this separation is equivalent to a 
high intensity of the driving laser field. This assump- 
tion allows for the parametrization of the peak only by 
the height and the width at half height. To calculate 
the peak height, we expand Eq. ( jl2| ) with respect to 
73/7 *C 1- The zcroth order corresponds to the usual 
Mollow spectrum without the additional peak as shown 
in Fig. |3| The first order corresponds to the additional 
peak. Thus the peak height A is given by the difference 
of the full spectrum and the Mollow spectrum at w = lul ■ 
We find 



A 



2(A 2 + 7 2 )ft 2 



(n 2 + 2a(2A 2 + 2 7 2 + Q 2 )) 2 7T 73 



(13) 



as an approximation for the peak amplitude A which is 
valid if the condition in Eq. ([!]) is satisfied. 

The peak width can be calculated without further 
approximation from the real part of the corresponding 
eigenvalues of the time evolution matrix B of the system 
PH , but the resulting expression is too complicated to 
be shown here. It turns out to be difficult to simplify the 
expressions for the width of the additional peak with pa- 
rameters chosen at the maximum peak intensity, where 
57, 7 and A are all approximately the same. It is however 
possible to make the Rabi frequency fl very large without 
loosing much of the peak intensity, if one appropriately 
changes the detuning A with the laser intensity. This can 
easily be done in experiments. This approach is similar 
to the secular limit usually made in explaining the Mol- 
low spectrum in the dressed state picture. However the 
elastic component and the additional peak are not secu- 
lar features and vanish for increasing Rabi frequencies if 
the other parameters are kept fixed. The reason for this 
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is that the driven atomic transition saturates such that 
the atomic dipole the laser field couples to vanishes. An 
adequate choice of A effectively prevents this saturation. 
We choose A such that the peak amplitude on the right 
hand side (rhs.) of Eq. (Oh has a maximum value: 



A 2 = 

max 



(1 + 2a)Q 2 - 8a 7 2 



(14) 



This is the analytical expression for the detuning A max 
for which we already obtaine d a nu merical value by read- 
ing off from Fig. ^in section III A. Our choice is justified 



by the fact that the variation of the amplitude in the pa- 
rameter range around the maximum intensity values are 
larger than the variations of the width. The absolute 
peak intensity is plotted against the Rabi-frequency and 
the detuning in Fig. (?]. For any fixed value of A, the 
intensity vanishes for large Rabi frequencies O. The con- 
dition Eq. ( |l4| ) is plotted as the solid lines on top of the 
saddle in the picture. It determines the detuning A such 
that for any given value of f2, the absolute peak intensity 
is close to its maximum value. In the approximation 



fl > 7, A = A, 



7 > 73,72 



as explained above we find 



pSCC 

peak 



A SCC 

peak 

rsec 
peak 



"73 



6a(l + 2a) 

l + 6a 
1 



9a(l + 2a) 7T73 

pscc a sec 
' peak ' peak 



3(1 + 6a) 



(15) 

(16) 

(17) 
(18) 
(19) 



for the width r p c c c ak , the amplitude ^4p C G c ak and the inten- 
sity /peak °f the narrow peak. Thus in this limit, the 
intensity of the additional peak only depends on a. A 
change of 72 and 73 with constant a only alters the shape 
of the peak, whereas the other parameters do not have 
any effect. These results agree with and explain the re- 
sults obtained by the numerical analysis. 

With the parameters in Eqn. ( |l6| ) - (|l9|), the contri- 
bution S pea k(ijj) of the narrow peak to the normalized 
fluorescence spectrum is 



Ipeak 



rscc 
peak 



pscc 
peak 



( w - to L f + (r p - k ) 2 



(20) 



Both the exact spectrum and a spectrum obtained using 
the above approximations are shown in Fig. |^. A closeup 
on the frequency range around u> = ujl is plotted in Fig. 
^. In both pictures, the solid line is given by the exact 
expression for the normalized inelastic spectrum in Eq. 
( Of ) . The dotted lines were calculated by adding the rhs. 
of Eq. (20) to the two-level Mollow spectrum obtained as 
explained above. It can be seen that the approximate pa- 
rameters in Eqn. (|l^) - ( |l9| ) are in very good agreement 
with the exact results. 



2. Random telegraph approach 

This subsection differs from the other parts of this pa- 
per in that it does not try to argue that the random tele- 
graph model is a good description of the physics of the 
additional narrow peak. Instead, it uses the model as the 
starting point for further calculations assuming that it is 
applicable. This will be justified in section IV A. Nev- 



ertheless, the fact that the results we will derive agree 
with other results in the respective limits gives further 
support for the model. 

Starting from the random telegraph model, the inten- 
sity of the additional peak can be calculated in a different 
way. According to the model, the additional peak stems 
from the stochastic modulation of the clastic component 
in the spectrum due to long dark and light periods in the 
resonance fluorescence signal. These modulations par- 
tially broaden the elastic component, but do not affect 
the relative weights of the inelastic Mollow components. 
The modulations only shift relative intensity from the 
elastic peak to the additional narrow peak. Therefore 
the sum of the relative intensity of the elastic component 
and of the additional peak has to be the same as the 
relative intensity of the elastic component in a system 
without these modulations. Thus the relative intensity 
of the additional peak I^ e ^ is equal to the difference of 
the relative intensities of the elastic component in a two- 
level system 1^1^ and the corresponding value i^stk; 1 
of the system discussed here: 



rphys r2-level 

peak ol 



clastic 



T-3-lcvcl 
clastic 



(21) 



The intensities of the elastic components are easy to cal- 
culate as in Eq. ( |Io| ) using the steady state values in Eq. 
(^) of the system and the corresponding values for the 
two-level system. Under the assumption that the ran- 
dom telegraph model is applicable, this yields as an ex- 
act expression for the relative intensity of the additional 
peak 

/P hys = 2(A 2 + 7 2 -4a 7 7 3 - 2a73 2 )r; 2 
P cak K{Q 2 + 2a(K + 4773 + 2 7 2 )) 1 ' 

where K = (2 A 2 + 27 2 + fi 2 ). This result is valid for any 
possible choice of parameters. 

If Eq. (0) is satisfied, a series expansion with respect 
to 73/7 up to zeroth order yields as an approximation for 
the peak intensity /p ea ^: 



rphys 



2(A 2 



m 2 



pcak ~ K(n 2 + 2aK) 



(23) 



Then the expression in Eq. (|T^) for the peak amplitude 
can be used to calculate an approximate expression for 
the peak width rp"^: 



^phys 
pcak 



f} 2 



2a(2A 2 + 2 7 2 + Q 2 ) 



2A 2 + 2 7 2 + Q 2 



■ 73 



(24) 
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Thus Eqn. (|l|), @ and (§4|) describe the amplitude, 
the intensity and the width of the peak if Eq. (Q) is 
fulfilled. In the limit in Eq. fll5| ), these results simplify 
to the previous results Eqn. (|16| ) - (|l9|). 

This method of calculating an exact expression for the 
peak intensity is not restricted to the atomic system dis- 
cussed here. It can be used in any system where the 
stochastic modulations due to light and dark states ac- 
count for a partial broadening of the elastic peak. For 
example, it is easy to derive an exact expression for the 
peak intensity in the atomic system driven by two lasers 
only by subtracting two steady state values using this 
technique. In the respective limit, this expression simpli- 
fies to the approximate result given in Q, Eq. (41). 



IV. PHYSICAL ORIGIN OF THE ADDITIONAL 
PEAK 

It is of course of great interest why there is an addi- 
tional structure in the resonance fluorescence spectrum 
compared to a two-level system. In several previous pub- 
lications, a random telegraph model was used to explain 
the dynamics in atomic systems like the one discussed 
here @ |l| ^ [l], g @ 0, j§. This 

model is applied to the present atomic system in the first 
subsection IV A . In the second subsection [V B , a spring 



model is used to give an intuitive picture of the dynamics 
of the atomic system. 



A. Random telegraph model 

If the fluorescence light of a system exhibiting the ad- 
ditional narrow peak is recorded by a broadband photo- 
detector, the signal shows a random sequence of bright 
and dark periods. Bright periods consist of a closely 
spaced sequence of fluorescence pulses, whereas dark 
phases are long periods in which no photons are emit- 
ted. The bright phases correspond to the usual atomic 
dynamics between the two levels connected by the driv- 
ing laser field as it is also present in a two-level system. 
The dark phases correspond to a trapping of the atomic 
population in a metastable state which in general is a 
coherent superposition of the atomic states. Because of 
the abrupt changes between light and dark phases, the 
fluorescence light behaves like a random telegraph signal. 

In this random telegraph model was used to ex- 
plain the occurrence of an analogous additional peak in a 
different atomic system. In their system, both the tran- 
sitions 1 <-> 3 and 1 <-> 2 are coherently driven by lasers, 
but no transitions from level 3 to 2 are allowed (73 = 0). 
The laser driving the transition 1 <-» 2 is assumed to be 
weak compared to the other laser, and 72 is set to zero 
during the calculations. Thus in the system in fl2|| , the 
additional atomic energy level |2) is coherently coupled 
to one of the other levels, whereas in our system it is in- 
coherently coupled to both driven levels by radiative de- 



cays. The result in jl^] was that the random fluctuations 
of the fluorescence light due to the change between light 
and dark states partially broaden the coherent spectral 
component. This broadening is essentially determined by 
the average length of the bright and dark periods. 

The fact that only a part of the elastic component be- 
comes broadened can easier be understood if one thinks 
of an ensemble rather than of a single atom. In fact, a dif- 
ferent approach to calculate the spectra has to be taken 
to describe single atom spectra |23|, |24|. In an ensemble, 
at any given time there is always a certain proportion of 
atoms in a light period. Because all the coherently emit- 
ted light is in phase with the incident laser, the sum of 
all these contributions forms the elastic part of the total 
fluorescence. This can be seen using the analytical ex- 
pressions for the average dark period length td and the 
average bright period length tb which will be derived in 



section 



©, (f§Mll) and (| 



IV A l|. If Eq. (|1| ) is satisfied, we have from Eqn. 



TB 



*coh 



rphys 
peak 



*coh 



TD 



7-phys 
peak 



TD + T B 



rphys 
peak 



*coh 



where (Ipe^k + ^coh) is the relative intensity of the co- 
herent peak in an unbroadened system not exhibiting 
the additional narrow peak. This intensity is equal to 
the relative elastic peak intensity in a two-level system. 
Thus a single atom is the proportion a of the time in a 
light period. Because of the phase coupling of all coher- 
ently emitted light to the driving laser field, in a large 
ensemble, this proportion of the total number of atoms 
can be considered as radiating continuously. Therefore 
effectively the proportion a of the unbroadened elastic 
intensity remains in the elastic part of a partially broad- 
ened spectrum. The remaining part f3 of the ensemble 
effectively can be seen as never having a bright period, 
but only dark periods. Thus the proportion (3 of the 
unbroadened elastic intensity is shifted to the additional 
peak. 

There has also been an experiment measuring the flu- 
orescence spectrum of single Yb(171) and Yb(172) ions 
]l4j . The spectrum of the first isotope exhibits dark and 
light periods, whereas the second has a continuous spec- 
trum. The atomic system used in this experiment is not 
the system presented in [Q, but the system described 
here with a fourth ancilla level (see Fig. ||). Still, the 
experimental results were reported to be consistent with 
the theory in Q . 

The mechanism for the exhibition of light and dark 
periods is not the same in the two atomic systems. In 
the system driven by two lasers, population coherently 
moves between the fluorescent transition and the addi- 
tional atomic state. Thus a dark state is not character- 
ized by all the population being trapped in the additional 
atomic level, but by a stable superposition of all three 
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atomic levels. This can be seen by the fact that a dark 
phase can end by an emission from both the upper and 
the metastable state fTsj] . In the system presented here, 
all the population is trapped in the additional state dur- 
ing a dark period; dark periods can only end by a transi- 
tion from this level |2(J. Because there are two different 
mechanisms which lead to analogous observations, it is 
reasonable that it is only the common fact that there are 
light and dark periods which accounts for the additional 
structure in the spectrum. Therefore the random tele- 
graph model should be a good theoretical model when- 
ever a system exhibits long bright and dark periods, no 
matter what the specific atomic structure is [ fL7| . 

As it seems, the experimental results were compared 
to results derived from the random telegraph model. The 
results of this model only depend on the average length 
of dark and light periods, Td and Tg. These parameters 
were not calculated using the experimental parameters, 
but measured during the experiment. However when td 
and tb are measured, one does not take into account 
the specific structure of the given atomic system, but 
only tests the random telegraph model. On this level of 
description, both the atomic systems with coherent and 
with incoherent coupling of the third atomic level are 
equivalent. This explains the results in As the mea- 
sured spectra are in quantitative agreement with the pre- 
dictions from the telegraph model, the experiment justi- 
fies the usage of this model. 

One way of testing the theoretical description of a spe- 
cific atomic system is to calculate the spectrum without 
using the telegraph model and to compare it to the mea- 
sured spectra. Another way is to calculate td and tb 
using a quantum jump approach as in Then the 

calculated durations can be compared to the measured 



ones. This will be done in paragraph IV A 1 



There is another argument for the random telegraph 
model. If this model is correct, then the fluctuations in 
the spectrum should not change the relative intensity of 
the Mollow spectrum, but only shift intensity from the 
elastic peak to the new narrow peak. A comparison of 
the relative intensity of the coherent part of a two-level 
system and of the sum of coherent and peak intensity 
in the system discussed here shows that this is true up 
to small errors due the approximations made to derive 
the intensity of the additional peak. Using this physical 
interpretation for the origin of the additional peak, one 
can calculate the peak intensity by taking the difference 
of the two-level coherent intensity and the three-level co- 
herent intensity, which can both b e calculated without 
approximations, (see section III B 2 ) 



The random telegraph model also allows for a dressed 
state description of the additional peak. In |^l| the oc- 
currence of long dark and light periods is explained in 
a dressed state picture. As the additional narrow peak 
stems from these long average durations, this description 
can be seen as a valid explanation for the narrow peak. 



1. Calculation of td and tb 

To calculate the average durations of bri ght and dark 
periods, we follow a method proposed in (g2|. At the 
time t = 0, we assume the atom not to be in the addi- 
tional atomic level, i.e. £22(0) = 0. We need to find the 
expectation value tb for the time the atom stays in the 
driven subspace {|1), |3)} without moving into |2). If at a 
later time t\ the atom is still in the driven subspace, the 
life expectancy at that time is again tb ■ If -P31 (ii ) is the 
probability of remaining in the driven subspace without 
moving into |2) during the time t\, we have 



T B = P 3 l(tl)(tl + T B ) + (1 - P 3 l(*l)) fh 



(25) 



where < / < 1. The first term represents the case 
in which the atom remains in the driven subspace. The 
second term is the probability of not staying in that sub- 
space multiplied by the time ft\ at which the atom moves 
out of the driven subspace. Thus we have 



TB = 



h(l-f) 



(26) 



For large t%, the atom may jump out of the subspace and 
back into it. This stops the current bright period, but 
does not change p22- For ti — *■ however the probability 
of leaving the subspace is dominated by single jumps into 
the additional atomic level. Therefore we have 



lim(l-P 3 i(ti))= limpaafr) (27) 

11 — >U £1 — >U 



and we obtain 



and thus 



tb = hm — 

tl-»0 p 2 2(tl) 



-1 ( dp22 



(28) 



(29) 



with the condition />22(0) = 0. A similar argumentation 
for the time the atom continuously stays in |2) yields 



'dp 2 



\ I "-P22 



(30) 



with p22{t = 0) = 1. 

From the equations of motion Eqn. (|J)-(||) we have for 

P22'- 

P22 = 273/O33 - 272P22 

This yields for the average length of a dark period using 
Eq. © 



-(273P33(0)- 2 72 p 2 2(0)) 



272 



(31) 



where we have used P22 (t = 0) = 1 and thus ^33 (t = 0) = 
0. This result is what one would expect, as the system is 
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FIG. 12: Relation between the bare and dressed states in the 
system of interest. The figure is plotted for the case A > 0. 
The solid arrows show the bare state transition widths; the 
dashed arrows denote energy splittings. (N), (N — 1) and 
(iV + 1) are the respective photon numbers in the laser field 
mode. f2 = \/Q 2 + A 2 is the generalized Rabi frequency. 



in the additional atomic level during a dark phase. Thus 
the length of such a dark period is given by the inverse 
of the total decay width 272 of this level. 

For the average length of a bright period we obtain 
from Eq. ( p9j ) under the condition p22(t = 0) = 0: 

t b x = 273P33(0)- 2 72jO22 (0) 
= 273,033(0) 

As the atom only evolves between the driven atomic lev- 
els, and as we look at the atomic system long after the 
laser was switched on when it has reached its quasista- 
tionary state, we can replace /933(D) by the steady state 
value for ^33 in the case of a two-level system. Using 
a = 72/73 and setting 73 to zero in Eq. (Q) yields: 



ss, 2 — level 
P33 



ft 2 



2(2A 2 + 2 7 2 + ft 2 ) 



Thus the result is 



ft 2 



2(2A 2 + 2 7 2 + ft 2 ) 



273 



(32) 



As one would expect, the average length of a bright pe- 
riod is inversely proportional to the total decay width 
2 7 3 out of the laser driven subspace. 

Equations Eqn. (|3l|) and (^) show that our general 
assumption in Eq. (fij) of a weak coupling of the addi- 
tional atomic level to the driven levels is equivalent to 
long average dark and bright phases. 

As shown in , one can derive an expression for the 
width r[ > e e 1 gg raph of the additional narrow peak in terms 
of the average length of dark and bright periods if these 



durations are long: 

ppeak _ -1 , -1 
1 telegraph B ' D 

Inserting the rhs. of Eqn. ([u]) and (|32| ) we get 
_ ft 2 + 2a(2A 2 + 2 7 2 + ft 2 ) 



^peak 
telegraph 



2A 2 + 2 7 2 



ft 2 



73 



(33) 



which is the result we already obtained in Eq. (g4|) . This 
is not surprising as Eq. (J2J) was derived under the as- 
sumption Eq. ([j]) which guarantees long average lengths 
of the dark and bright periods. 

If the parameters describing the atomic system and the 
laser field are known, expressions Eqn. ( |3l| ) and (|3^) can 
be compared to average light and dark period times ob- 
tained by experiments as Q. Alternatively, Eq. ( |33| ) 
may be compared with an experimentally determined 
width of the narrow feature. This allows to test the 
quantum optical description of the atomic system via the 
telegraph model. Equally assuming that this description 
is appropriate, these results can for example be used to 
measure the decay widths 72, 73 or the effective decay 
width from |2) to |1) of a system of ancilla pumping laser 
and decay from |4) to |1) as in Fig. 0. 



B. Spring model 

In this section, we will use the picture of a simple spring 
model to describe the atomic dynamics. To physically in- 
terpret the calculated spectrum, we transfer the system 
to the dressed state picture (c.f. Fig. |l^). The dressed 
states are defined as the eigenstates of the system Hamil- 
tonian in the interaction picture in Eq. (E2): 



with 



N) 
N) 



sin(0) 
cos(0) 



sin(0) |3, N - 1) + cos(0) |1, N) 
cos(0) |3, N-l) -sin(0) |1, N) 



ft 

^(A 2 + ft 2 - AVA 2 + d 2 ")! 



1 -sin 2 (0) 



(34) 



(35) 



N and N — 1 are the respective photon numbers in the 
laser field mode. The third bare state |2) is not coupled 
to the driving laser field and thus remains the same in 
the dressed state picture. 

By choosing a detuning according to Eq. (14), the 
Rabi frequency of the driving laser field can be made very 
large as compared to all transition widths in the system 
without loosing weight in the elastic part of the spectrum 
or the additional peak (see Fig. |?|). This secular limit 
allows for a separation of the central components whose 
time evolution is at the laser frequency u>l and of the 
sideband components. 
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FIG. 13: The spring system used as a model for the atomic 
dynamics. 



The resulting equations for the corresponding popu- 
lation vector g of the central frequency parts at ojl in 
secular approximation can be written as 



d - ™ - 
- e = M Q 



(36) 



[Q++, Q — , Q22) 



Qu for i e {+, — , 2} are the populations of the dressed 
states summed over all possible laser field photon num- 
bers. 
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FIG. 14: Fundamental modes of the spring model, (a) corre- 
sponds to the elastic component, (b) to the Mollow peak, and 
(c) to the additional narrow peak. The arrows show the rel- 
ative elongation directions of the masses. The small arrow in 
mode (b) is a correction term only used in the closer analysis. 



The time evolution matrix M is given by 



M 



-2 73 sin 2 0- 2 7 sin 4 
2 7 sin 4 
273 sin 2 



27cos 4 272 cos 2 

-273 cos 2 - 27 cos 4 272 sin 2 
273 cos 2 —272 



An essentially equivalent picture for these spectral 
components can be obtained using a mechanical spring 
model as shown in Fig. 



this spring model are 



13. The equations of motion for 



dt 2 



Q = F-Q 

Q = {m + q + ,m-q-,m 2 q2) T 



(37) 
(38) 



with 



-D1+-D2 

-Pi 

771-1- 
771-1- 



Dl±D3 

rri — 
D 3 



D 3 
m 2 
D2+D3 



rrii are the masses and (i € {+,—,2}) are the elon- 
gations of the masses. Di (i £ {1,2,3}) are the spring 
constants as denned in Fig . [T3| . 

After transforming Eq. (|36D to a second order differen- 
tial equation in t and identifying the atomic population 
with the product of mass and elongation, a comparison 
with the spring model equation Eq. (|37|) yields the rela- 
tions 



m 2 



sin- 4 (0) 
cos - 

73 

72 



\0) 

sin- 2 (6>) cos -2 (0) 



Do 



37 2 + 4773 - 47273 + 7 2 cos(40) 

73 

■ 2 {7 + 4 72 + 273 
8X0.(0) 

-2(7 + 73) cos(26») + 7 cos(46i)} 

73 



-{7 + 472 + 273 



cos 2 (0) 

+2(7 + 73) cos(26») + 7 cos(46»)} 

for the masses and the spring constants. 

From the above equation for mi we can see that the 
parameter a — 72/73 determines 7712 which is the mass 
corresponding to the additional atomic level. This ex- 
plains why a is a crucial parameter for the system. For 
example in the secular limit, the peak properties depend 
on a although all other parameters have only very little 
influence on the results as shown in the analytical calcu- 
lations. 

There are three fundamental oscillation modes in the 
spring model, which can be calculated from the eigen- 
vectors of the time evolution matrix F. See Fig. [l4| 
for a symbolic depiction of the modes. The arrows repre- 
sent the relative elongation directions of the masses. Fig. 
( |l5| ) shows the possible transitions between the dressed 
states with the respective transition widths. We denote 
the transition rate from i to j as € {+, — , 2}). 

The first mode (a) corresponds to a uniform motion 
of all three masses. As there is no relative movement 
between the masses, none of the spring constants - or 
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FIG. 15: Possible transitions between the dressed states in 
the spring model with corresponding rates. 



in the atomic model none of the decay constants - are 
relevant. Thus the total decay width r| pring in this mode 
is zero: 

pspring = q 

This mode corresponds to the elastic peak which is of 
zero width. 

In the second mode (b), the masses + and — oscillate 
against each other, whereas mass 2 is approximately at 
rest. We begin by neglecting the small movement of mass 
2 shown in Fig. which will be included later in a 
closer analysis. In this approximation, only the spring 
constants or decays between + and — are important for 
the time evolution of the system, which are of order of 7 
(c.f. Fig. [l5|). Therefore this mode corresponds to the 
inelastic Mollow peak at oj = ujl. The expected decay 
width r b pimg is the sum of the decay widths from + to 
— and vice versa: 



-psprmg 
1 b 



= 2 7 sin 4 (0) + 2 7 cos 4 (0) 

3 + cos(4(9) 
= n 7 



The exact width Tb of this peak calculated in secular ap- 
proximation from the equations of motion for our three- 
level system Eq. ( |36| ) up to first order in 73/7 is 

3 + cos(40) 5 + 3cos(40) 
Fb = ^^ + 3 + cos(W 73 ^ 

Thus the predicted and the exact widths are the same up 
to a correction term of first order in 73/7. 

In the third mode (c), + and — oscillate against 2. 
Only the transitions from + and — to 2 and vice versa 
are used, but not the ones between + and — . As the 
transitions between 2 and + or — are of order 73, this 
mode corresponds to the additional peak. The expected 
width r^ rin s is 

= r+2 + r_ 2 + r 2+ + r 2 _ 

= 2 73 sin 2 (0) + 2 73 cos 2 (6>) 

+27 2 cos 2 (0) + 2 72 sin 2 (0) 
= 2(a+l) 73 



The exact width T c of this peak in secular approximation 
from Eq. ( |36| ) up to first order in 73/7 is 



r c = 2(a + l) 73 - 



5 + 3cos(40) 
3 + cos(40) 



73 



(40) 



Also for this mode the widths are the same up to a cor- 
rection term of first order in 73/7. 

Using this simple picture, it is possible to predict the 
widths of the spectral components up to a correction term 
of order 73 for the last two modes. This is a minor cor- 
rection for the Mollow peak, which is of order 7, but can 
be a large correction to the width of the additional peak, 
which is also of order 73. It is interesting to note that 
the correction term in both modes is the same up to a 
change in sign. 

A closer analysis shows that the populations of the 
dressed states in the modes corresponding to the Mollow 
and the additional peak are not the same. Therefore 
the relevant transition widths do not contribute equally 
to the total width of the respective modes. This is the 
reason for the correction terms found above. Thus in the 
following, we use the populations of the dressed states 
in the different fundamental modes as weights for the 
transition widths and take the small correction due to 
the movement of mass 2 in the Mollow peak mode into 
account. By this it is possible to obtain the exact secular 
approximation results Tb and T c in first order in 73 /j for 
the widths of the Mollow and the additional peak from 
the spring model. 

For mode (b) corresponding to the central inelastic 
Mollow peak we have as relative weights in the dressed 
states up to first order in 73/7 from the corresponding 
eigenvector: 





a h 


w b = 


2 




«b 


w b = 


2 




Ob 


w l = 


2 



= ^ 1- 



m_ — m+ 



m+ + m_ 
- m + 73 
m + + rri- 7 



J3_ 

1 



in 



"b" denotes the mode, and the super- 
and "2" mark the corresponding masses 



The subindex 
scripts "+", "- 
or dressed states. a\, — + + w 2 is a normalization 
constant. This constant is due to the fact that the eigen- 
vector from which the weights are calculated does not 
have a definite norm, such that only the ratios between 
the different weights can be calculated from it. The nor- 
malization constant does not have to be unity, as only 
the total population has to be normalized to unity. This 
total population is the solution to Eq. (|3^) and con- 
sists of a linear combination of the contributions from 
the three different modes, thus allowing for a different 
normalization of the weights of a single mode. 

Clearly the weights should not be larger than unity, 
because the transition widths cannot contribute with a 
prefactor larger than one to the total width of the peak. 
Thus we have < ab < 2. As in a peak due to a single 
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decay the corresponding single transition width is equal 
to the total width and as additional transitions can only 
increase the total width, it is plausible that the largest 
weight is unity. Therefore we assume ab = 2. 

The relative weights of the dressed states in mode (c) 
corresponding to the additional peak up to zeroth order 
in 73/7 are 

+ u c m+ 

wZ = — ■ 

2 m + + m_ 

a c m_ 
w c = — ■ 

2 m + + ?n_ 

2 Q c 

w r = — 
2 

Again, the subindex "c" denotes the mode. A similar 
argument as for mode (b) yields a c — 2 as assumption 
for the normalization constant. It is not necessary to take 
higher orders in 73/7 into account for the narrow peak 
mode weights as the corresponding transition widths for 
the mode (c) are already of first order in 73/7. 

With these expressions for the weights, we obtain as 
width of the Mollow peak up to first order in 73/7: 



w b ( r +- + r +2) + < r_ 

3 + cos(46>) _ t 5 + 3cos(4<9) 
3 + cos(46>) 



■7- 



73 



The width of the narrow peak in first order of 73/7 
calculates to 



+2 



2(a+l)73 



w c r_ 2 + w c (r 2+ + r 2 _ 

5 + 3 cos(40) 



3 + cos(46») 



73 



(41) 



These results agree with the exact secular approximation 
results Tb and T c in first order in 73/7, Eqn. ( |3^ ) and 
©. This justifies our choice for ab and a c . 

In the secular limit, the result for the width of the ad- 
ditional peak is also in agreement with the corresponding 
expressions obtained in the analytical section, Eqn. (24) 
and (33). Inserting the explicit values Eqn. (|34|) and (3£) 
for the trigonometric functions into the rhs. of Eq. (41) 
yields 



pin 



n 2 + 2a(2A 2 + n 2 

2A 2 + n 2 



73 



Taking the secular limit in the expression for the peak 
width Eqn. ( gjj ) and ( p3| ) calculated in the analytical 
section by expanding to zeroth order in 7/fi yields the 
same result. Thus the spring model is an adequate way 
of visualizing the atomic dynamics responsible for the 
central spectral features in the secular approximation. 

Of course also the satellite peaks in the Mollow spec- 
trum can be discussed in this way. But as there is always 
only one spectral component at each side, this does not 
give new insights. 

A corresponding spring model yields the correct re- 
sults also for a two-level system driven by one laser field. 



Thus the model is not limited to the three-level system 
discussed here and could be useful in predicting the pos- 
sible modes in more complex atomic systems. 

One might also think of extending the mechanical 
spring system in a way which allows to drive the oscilla- 
tions. As discussed in almost any textbook on mechan- 
ics, a suitably chosen driving frequency may give rise to 
resonance effects. Under ideal conditions, it is even pos- 
sible to excite only one of the modes of the spring model. 
As each mode corresponds to one of the spectral compo- 
nents in the atomic system, one might hope to eliminate 
all spectral components but one from the spectrum under 
these ideal conditions. There are two different ways to 
find a suitable spring model. On the one hand, one can 
try to modify the atomic system, find the corresponding 
spring model, and then check, whether the modifications 
allow to drive the spring system. On the other hand, 
one can alter the spring system such that it allows for an 
external driving and then try to translate it back to an 
atomic system. However a first attempt to modify the 
atomic system by an additional laser field on one of the 
transitions did not succeed to allow for a driving of the 
spring system. There are also some general arguments 
against the possibility to excite the various spectral com- 
ponents using this method. The spring model does not 
allow one to easily obtain the amplitudes of the respective 
modes. Therefore it is not clear whether a large excita- 
tion of a spring mode corresponds to a large intensity of 
the corresponding spectral component. Also, the masses 
and spring constants of the spring model depend on the 
parameters of the atomic system, and therefore also the 
resonance frequencies. Thus the parameters of the spring 
system corresponding to an atomic system continuously 
change during the excitation process, which is different 
to conventional spring systems. Nevertheless finding an 
atomic system which translates to a driven spring system 
may provide novel insights in the atomic dynamics. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the resonance fluo- 
rescence spectrum of a three - level atomic system driven 
by a single monochromatic laser field. The third atomic 
level adds an additional decay path from the upper laser 
driven state to the lower one by radiative decays. An 
additional narrow peak compared to the usual two-level 
fluorescence spectrum was found. By a numerical analy- 
sis, the peak was found to be especially pronounced for a 
very small coupling of the additional atomic level to the 
laser driven ones. Also the Rabi frequency of the driving 
laser must not be too high if the detuning of the driving 
laser is kept fixed. 

These results were used to derive analytical expres- 
sions for the width, the height and the intensity of the 
additional peak using two different methods. The first 
method uses a secular approximation in which a suitable 
choice of the detuning allows to increase the Rabi fre- 
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quency without losing intensity in the additional peak. 
In this limit, the peak intensity essentially only depends 
on the relation a — 72/73 between the decay width from 
and to the additional level, but not on other variables 
such as the Rabi frequency ft or the absolute value of 
the decay widths. In addition to the dependence on a, 
the peak amplitude is explicitly proportional to the decay 
width 73, and the peak width is proportional to 7 3 ~ 1 . In 
the second method, results from the following discussion 
of the physical origins of the peak were used to derive 
an exact expression for the intensity of the narrow peak 
under the assumption that the random telegraph model 
is appropriate to describe it. From this expression for 
the intensity, the width and the amplitude of the peak 
were calculated under the assumption of a weak coupling 
of the additional atomic level to the driven ones. In the 
respective limits, the results of the two different methods 
are equal. 

In the second part of the paper, we have discussed the 
physical origin of the additional narrow peak. We have 
applied the random telegraph model to our atomic sys- 
tem. In previous works, this model was used to explain 
the exhibition of similar additional spectral features in a 
different atomic system. In this alternative system, the 
three atomic levels are coupled by two lasers in a V - 
configuration instead of the additional decay path in our 
system. The average length of light and dark periods, 
which are the crucial variables in this model, were calcu- 



lated analytically. Using these results, we have explained 
the relation between the atomic system discussed here, 
the atomic system where both transitions are driven by 
laser fields, and an experiment which was recently re- 
ported. In the experiment, a setup as described in this 
paper was used. The experimental results are in quan- 
titative agreement with the theoretical predictions from 
the random telegraph model. 

Finally we have discussed a simple spring model to 
explain the atomic dynamics in the secular limit. This 
model provides an intuitive picture for the different spec- 
tral contributions in the fluorescence spectrum and allows 
to visualize the atomic dynamics responsible for the flu- 
orescence spectra. Also it can be used to easily calculate 
approximate expressions for the widths of the different 
peaks. 
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